Multipartite entanglement in fermionic systems via a geometric 

measure 

Behzad Lari*, P. Durganandini^ and Pramod S. Joag* 
Department of Physics, University of Pune, Pune Jj.ll 007, India. 

Abstract 

We study multipartite entanglement in a system consisting of indistinguishable fermions. 
Specifically, we have proposed a geometric entanglement measure for N spin-^ fermions dis- 
tributed over 2L modes (single particle states). The measure is defined on the 2L qubit space 
isomorphic to the Fock space for 2L single particle states. This entanglement measure is de- 

W ■ fined for a given partition of 2L modes containing m > 2 subsets. Thus this measure applies 

00 . . 

to m < 2L partite fermionic system where L is any finite number, giving the number of sites. 

i-^ | The Hilbert spaces associated with these subsets may have different dimensions. Further, we 
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have defined the local quantum operations with respect to a given partition of modes. This 
definition is generic and unifies different ways of dividing a fermionic system into subsystems. 
We have shown, using a representative case, that the geometric measure is invariant under local 



to calculate multipartite entanglement in some correlated electron systems. To the best of our 



unitaries corresponding to a given partition. We explicitly demonstrate the use of the measure 
00 

O 

knowledge, there is no usable entanglement measure of m > 3 partite fermionic systems in the 
literature, so that this is the first measure of multipartite entanglement for fermionic systems 



going beyond the bipartite and tripartite cases. 
PACS numbers: 03.65.Ud, 71.10.Fd, 03.67Mn, 05.30Fk 
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INTRODUCTION 



The nonlocal correlations implied by the states of multipartite quantum systems form 
the basis of quantum information processing technologies like quantum teleportation and 
quantum computation [lH3(. Recently the question of understanding and using entan- 
glement in the systems of identical and indistinguishable particles has seen a surge of 
interest [4H7|. Two or more identical particles become indistinguishable when their wave 
functions overlap. Such a situation can arise, for example, in a quantum device based on 
quantum dot technology [8|, |9|. Here qubits are realized by the spins of the electrons in a 
system of quantum dots. The overlap between the electron wave functions in different dots 
can be varied by controlling parameters like gate voltages or magnetic field which change 
the tunneling amplitudes of the electrons from one dot to the other. For non-negligible 
overlaps, the entanglement between the qubits is then intimately connected to the elec- 
tron entanglement which is essentially that of indistinguishable fermions. Entanglement is 
also expected to play a fundamental role in many physical phenomena like quantum phase 
transitions, quantum Hall effect, etc 10|, |ll| involving many body quantum systems. 

The study of entanglement in many body quantum systems consisting of identical 
particles has posed challenging fundamental questions about the definition and nature of 
entanglement in such systems due to the inherent indistinguishability of particles 4,112- 
ll4J . For example, the (anti)symmetrization necessary for indistinguishable (fermions) 
bosons already leads to quantum entanglement. However, it is known that the correlations 
due to symmetrization or anti-symmetrization are not by themselves a physically useful 
resource for quantum information and communication technologies; for example, there 
is no measurement we can do locally on a fermion in a localized state which is affected 
by the existence of identical fermions in other parts of the universe 151 ] . Most studies 
of entanglement in many body systems consisting of identical particles have focussed 
on the study of bi par tite entanglement. There have been few studies on multipartite 



entanglement 



lfl 



17j | although the study of multipartite entanglement in such systems 



is also an important and interesting question. 



In this work, we study multipartite entanglement in a system consisting of identi- 
cal particles. We use the idea due to Zanardi |4j whereby the Fock space of a system 
of fermions is mapped to the isomorphic qubit or 'mode' space. We then discuss en- 
tanglement in this 'mode' space via a geometric measure. The idea is to use the Bloch 



representation of the state of the m-partite quantum system [181 ] . The measure is defined 
by the Euclidean norm of the m-partite correlation tensor in the Bloch representation. 
This correlation tensor contains all information of genuine m-partite entanglement (see 
section III). Such a measure was proposed earlier by one of the authors (and Hassan) for 



N qubit and N qudit pure states 19|, |20| and shown to satisfy most of the properties 



expected of a good measure. The Bloch representation of a quantum state has a natural 



geometric interpretation, which is why we call this measure a geometric measure 18] 



Other geometric measures are based on the distance of the given state from the set of 



separable states in the Hilbert space 
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22]. An important question in the context of 



quantum entanglement is that of locality. For indistinguishable particles distributed over 
'modes' (which are taken to be single particle states of particles constituting the system) , 
local operations have meaning only in the context of partitions over modes. A quantum 
operation confined to a single subset in a given partition is then a local operation. We 
therefore define entanglement in such a system with respect to partitions and require it to 
be invariant under local unitaries defined with respect to a given partition. We explicitly 
demonstrate the use of the measure to calculate entanglement in some correlated electron 
systems. 

The outline of the paper is as follows: In SecJTTl we begin by briefly reviewing some 
details about the isomorphism between the Fock space of a system of indistinguishable 
particles and the 'mode' space |4j. This will help us in defining various quantities and also 
set up notation necessary for the subsequent analysis. In Sec lHIl we define and construct 
the geometric measure. Various properties of the measure are also discussed with reference 
to a specific example in Sec JIVI In SecJVJ we use the measure to study entanglement in 
the Hubbard dimer and trimer. Finally, we conclude in Sec JVIl 



II. MAPPING BETWEEN FOCK SPACE AND QUDIT SPACE 

We deal with N spin-| fermions on a L site lattice. The total number of available 
(localized) single particle states are then 2L in number. The fermionic Fock space in 
the occupation number representation has basis states of the form \n\n 2 . . . n 2 i) (ni = 
0,1; i = l,..., 2L). We further assume that the total number of particles is conserved. 
This means that we only deal with subspaces of the Fock space corresponding to a fixed 
eigenvalue for the total number operator. We shall refer to this number super-selection 
rule as N-SSR. For a iV-fermion system, we call such a subspace of the Fock space l N- 
sector' and denote it by F/y. The iV-sector of a 2L mode system is the subspace Fn 
of dimension ( „) of the Fock space with the dimension of the Fock space for 2L single 
particle states (N = 0,1,2, ... , 2L) being 

Since a 2L qubit Hilbert space (C 2 )® 2L has exactly this dimension, it is possible to con- 
struct an isomorphism between the Fock space and the 2L qubit Hilbert space (C 2 )® 2L 



231 ] . The particular isomorphism we implement is 

|mn 2 . . . n 2L ) -> | ni) <g> \n 2 ) ® ■ ■ ■ ® \n 2L ) ; m = 0, 1 ; i = 1, 2, . . . , 2L (2) 

where, in qubit space we associate |0) -H- | f) and |1) -H- | l). Note that the Slater rank of 
the Fock basis states \nin 2 ■ ■ ■ n 2 L} is 1 so that these are separable states. Thus the above 
isomorphism maps separable basis states in Fock space to the separable basis states in 
qubit space. Further, the subspace structure of the Fock space namely, 
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F N (3) 



N=0 



is carried over to the qubit space under mapping (Eq.©) because each subspace of the 
Fock space with conserved fermion number iV is mapped to a subspace of the 2L qubit 
space spanned by the basis vectors with iV ones and 2L — N zeros. We can write 



1L 



H 2L = (C 2 f 2L = H 2L (N) (4) 



N=0 



where H 2L (N) is the image of F N in Eq.flH]) under the map given by Eq.(T2]). 

Next crucial step is to transfer the action of the creation and annihilation operators 
on Fock space to the qubit space, under the isomorphism given by Eq.(2) [24]]. We need 
the creation and annihilation operators a and a) acting on a single qubit state, 

o|0) = 0, o|l) = |0) 
a t|0) = |1), a f |l) =0 



(5) 



such that, 



a, ->■ I ® I ® ■ ■ - ® a <g> ■ • ■ <g> / 

ith qubit 
oj ->■ I ® • • • ® O^ g> • ■ • g> J 



(6) 



Here Oj (a]) is the annihilation (creation) operator acting on Fock space F 2 l, annihilating 
(creating) a fermion in ith mode. I is the identity on single qubit space. The tensor prod- 
uct satisfying the correspondence in Eq.([6]) must be consistent with the anti-commutation 
property of the Fock space creation and annihilation operators, 

{a h a)} = 6 i:j {a u a,-} = = {a], a]} (7) 

This requirement leads to the following action of the tensor product operators on the 2L 
qubit states 

(J®/®---® a(a r ) ®---®/)(|m}®---® |7ii)<8>---® \n 2L )) = (-l)£?=*+i n ' 

ith place 

(|ni) ® • ■ ■ <g> a{a))\rii) ® ■ ■ ■ ® \n 2L )) 

ith qubit 



Here rii e {0, 1} ; i G {1, 2, . . . , 2L} and J2j=i+i n j * s evaluated mod 2. Using Eq.©, it 



is straightforward to see that 

{J®/®- • -<g> a(a ] ) ® ■ • -®I , J®/®- • -(g) a(a f ) ® • • -®.f}(|ni)®- • -®|n;)®- • -®|rij)®- • -E^l)) = 

jt/i place jth place 



(9) 



and 



{J ® 7 ® ■ ■ • ® a, ®-- 

ith place 

(1%) ® • • • ® |nj) ® • • • ® |rij) 
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(10) 



We note that the phase 
of the parity operator 



actors appearing in Eq.([8]) are the consequence of the conservation 

m 

P = U^ 1 (l-2a\a l ). (11) 

Henceforth, in this paper, by 'fermions' we mean spin-| fermions. Further, we call a 
single particle state a mode. Thus two spin-^ fermions on two sites is a four mode system. 
In general, N spin-| fermions on L sites is equivalent to N fermions on K = 2L modes. For 
example, two spin-| fermions on a two site lattice, A, B say, generate four single particle 
states or modes \A ^), \A \), \B "\), \B \). In this work, we deal with entanglement between 
subsets forming a partition of a 2L-mode fermionic system. We define the entanglement 
measure for any such partition of a 2L-mode system without any restriction on the number 
and the size of the subsets forming the partition. These subsystems may involve different 
degrees of freedom, for example, we can deal with the entanglement between spins and 
sites or entanglement between two spins on the same site (intra-site entanglement). Or if 
each of the subsets partitioning the 2L modes comprises modes with common site label 
we have the entanglement between sites or the 'site entanglement'. Thus all physically 
realizable subsystems of iV fermions over 2L single particle states can be addressed by 
dividing the 2L modes into suitable partitions, for example the 'particle entanglement' 



is 



defined in |6|. 

We now define the local and non-local operations on the 2L mode fermionic system 
j#i6|, [26|, |27J. We do this by using the corresponding operations on the isomorphic qubit 



space H 2L . We note that, due to isomorphism between F 2L (Eq.(j2J)) and H 2L (Eq. (J5J)), 
partitioning 2L modes is equivalent to the corresponding partitioning of the 2L qubit 
system into subsystems. Locality is defined with respect to the partition of 2L qubits 
(or, the corresponding partition of 2L modes) between whose subsets we are seeking 
entanglement. The operations on the state space of a single subset in a partition of 2L 
qubits is taken to be local. The entanglement measure defined with respect to a partition 
must be invariant under a unitary operation which is local with respect to that partition. 
We will illustrate this point later, using the geometric entanglement measure defined 
below. Henceforth 'mode' and 'qubit' are taken to be synonymous and we shall use the 
expression 'modes' instead of 'qubits'. In other words, the spaces F 2 l and H 2 l are taken 
to be the same. 

III. GEOMETRIC MEASURE FOR ENTANGLEMENT. 

A. Definition 

We define a geometric measure for the partitions of the 2L mode N fermion systems in 
pure states. Although the definition of the measure is quite general, the fermion number 
super- select ion rule restricts the pure states to the appropriate subspace corresponding to 
N fermions, namely H 2 l(N) (Eq.(j4|)). Thus a state |0110) G H^{2) may be partitioned 
as 1 01) £g> 1 10) or as |011) <8> |0) etc where the definition of the tensor product is consistent 



with Eq. (El). We use 
to get this measure 



;he Bloch representation of N— partite states (drawn from H 2 l(N)) 
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20|. 



First we assume that a partition equally divides 2L modes into subsets, ie all subsets 
in the partition contain equal number of modes, say n. This corresponds to the case of 
H 2 l divided into subspaces of dimension d = 2 n , n being some divisor of the number 
of modes 2L. To get the entanglement measure we expand the state p = \ip){ip\ of the 
system, supported in the appropriate H 2 l(N), in its Bloch representation. 

In order to give the Bloch representation of a density operator acting on the Hilbert 



space C d <S> C d ® ■ • • <E> C d of an m = 2L/n-qudit quantum system, we introduce following 
notation |20]. We use k, ki (i — 1, 2, . . .) to denote a qudit chosen from m qudits, so that 
k, ki (2 = 1,2,.. .) take values in the set M = {1, 2, ... , m}. The variables a^ or a^ for a 
given k or /c; span the set of generators of SU(d) group for the kth or kith qudit, namely 
the set {Ai, A2, • • • , A d 2_ 1 } for the kith qudit. For two qudits k± and k 2 we define 



A£ = (Id ® h ® • • • ® A a ® J d ® • • ■ (8) I, 



^ = ^ ^ /d ® ' ' ' ^ A °*3 ® /d ® " " " ^ /rf ) 

^S^S? = ( J rf ® J d ® • • • ® A « fcl ® ^ ® • • • ® A ^ 2 ® J <* ® J d) 

where A Qfc and X ak occur at the fcith and k 2 th places (corresponding to fcith and k 2 th 
qudits respectively) in the tensor product and are the o^th and a^ 2 th generators of 
SU(d), a^ 12 = 1, 2, . . . , d 2 — 1 . Then we can write 

V V t \ (/ci) A (/C2) ---A (/cm) +---+ V i A (1) A (2) ---A (m) l 

{fci.fca,— .fcjvf} a fci a fe 2 '" afe Af aiaa— ow 

(12) 
where s( fc ) is a Bloch vector corresponding to fcth qudit, s( fc ) = [s a Jf ~\ which is a tensor 
of order one defined by 

s ak = ~Tr[p\<£] = d -Tr[p k \ a ^ 

where p k is the reduced density matrix for the kth qudit. Here {ki, k 2 , . . . , & A ,/}, 2 < M < 
m, is a subset of Af and can be chosen in (™) ways, contributing (™) terms in the sum 
EE fc k ... k > i n Eq. (1III A[) . each containing a tensor of order M. The total number of 
terms in the Bloch representation of p is 2 m . We denote the tensors occurring in the sum 
E {fcl>fc2 ,.., feM} , (2 < M < m) by r** 1 '**'"'*^ = [*«*, «*,-«*„] which are defined b ^ 

t = — TrfoA (fcl) A (fc2) ■ • • A (fcM) l 

L a kl a k2 ...a kM 2^ ^ a ki °*2 a *M J 



M 



where pk 1 k 2 ...k M * s the reduced density matrix for the subsystem {fci, fc 2 , . . . , Um}- Each 
of the (J^) tensors of order M, occurring in the Bloch representation of p, contains 
all information about entanglement of the corresponding set of M subsystems. All 
information on the entanglement contained in p is coded in the tensors occurring in 
the Bloch representation of p. The tensor in last term in Eq.( 1HI Al) . we call it T^ m \ 
contains all the information of genuine m-partite entanglement. This follows from the 
observation that all other terms in the Bloch representation of p (Eq. (1III Al) ) correspond 
to subsystems comprising M < m qudits and the density operator contains all possible 
information about the state of the system. 



The operators X ak , «& = 1, 2, . . . , d 2 — 1 are given by [29|] 



with 



A = {u 12 , u 13 , m 23 , . . . , v 12 , #13, v 23 , . . . , w t , t&i, . . . , u> rf _i} 

Vjk = —i(Pjk — Pkj) 



(13) 



wi = ^ I l(l + 1) (Pu + --- + Pu- IPi+i,i+i) 
!<]<k<d; \<l<d-l 



(14) 



where 



Ai=|*>(*| (k,l = l,2,...,d). 



Note that each of the generators of the SU(d) group Aj , i = 1, 2, . . . , d 2 — 1 acts on a 
single qudit space and hence is local (see SecHI]), apart from the phase factor contributed 



by their action, as given by Eq.flS]). We assume these phase factors to be absorbed in the 
coefficients in the expansion of the density operator p. 

Let a 1L mode N fermion system be partitioned by m = 2L/n subsets, each containing 
n modes. Then for this partition, we define the entanglement measure for a state \ip) G 
H 2L (N) by Id 



E 



i »<'p 



where 



\ a\—a m =\ 

and ||r|| se p is ||r|| for separable (product) m qudit state 



£ 



t 2 

"ax—an 



(15) 



(16) 



i *<'i< 



d{d-l] 



m/2 



(17) 



B. Entanglement in partitions with unequal subsets 



We can also generalize the definition of the entanglement to the case where the cor- 
responding qubit subsystems have unequal dimensions. We discuss the simplest case of 
bi-partite entanglement with partitions having unequal dimensions, say d\ and di- In this 
case, the definition of the geometric entanglement measure generalizes to 



E 



i »('p 



where 



\ 



d?-ld§-l 

E£«5 



where 



''./ 



I A,: ® A^ 



K, 



(18) 



(19) 



(20) 



Here \ (i — 1, . . . , d\ — 1) and Xj (J = 1, 
SU{d2) respectively. ||T sep || is given by 



d\ — 1) are the generators of SXJ{d\) and 



' sep\ 



,(1)1 



,(2) I 



dx{d\ - 1)\ fd 2 {d 2 - 1 



(21) 
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Here | \s^ \ | and | \s^ \ | are the norms of the Bloch vectors of the reduced density operators 
for each subsystem. 

It is straightforward to extend these definitions to partitions containing more than 
two subsets. 

IV. ENTANGLEMENT IN A FOUR MODE SYSTEM 

With the entanglement measure defined as above, we give an example wherein we can 
compute the entanglement for different partitions and illustrate our comments on local 
and non-local operations. We also compute the upper bounds on the entanglement. 

A. Local and Non-local operations 

Consider a four mode system and the normalized state \ip) G H^{2) defined as 

\ip) = -^{ia\1100} + |1001) + |0110) + |0011) + /3|0101) + |1010)} a 2 + (3 2 = 2 (22) 
V6 

where a, /3 are real. Note that \ip) can be treated as a member of the Fock space F^{2) 
with the kets appearing in it being its basis states. Consider the evolution of the system 
in state \ip) G F±(2) via the Hamiltonian 

H = /(0J04 + a^ai) + qh\fi2 + Thi + 7^3 + ^(a|a 2 + a\ai) (23) 

acting on F4. Here / term is the interaction between two modes on different sites (inter- 
site interaction), 77 term is the interaction between two modes on the same site (intra-site 
interaction). T and 7 correspond to single mode on site A and B respectively, q term 
involves number operators hi = ajai ; i — 1,2 for first two modes, on A site. We have 
included all the different kinds of typical interactions encountered in condensed matter 
systems, respecting number super-selection rule. After an infinitesimal unitary evolution 
via this Hamiltonian, the state 1^) evolves to 



M 



W) = \i/>)-ieH\1>) (24) 

11 



By employing the mapping of annihilation and creation operators in Eq.flf)]) and Eq.(jS]) 
and that of Fock space basis states in Eq.flSJ), we get, for \ip') 

IVO = A={(ia + ief + age) 1 1100) + (1 - ire - ier//3)|1001) + (1 - ije - ier/)|0110) + 
V6 

+(1 - ief - ze7)|0011) + (P - effi - ier])\0101) + (1 + ief - ieT - ie 7 - ier/)|1010) 

(25) 

Now we find the entanglement for different partitions of this four mode system, using 
the geometric entanglement measure, Eq. ffT5|) . We first partition four modes into four 
subsets, each containing one mode. This case gives genuine entanglement between four 
modes, which is more general than only the bipartite entanglement considered in the 
literature. For this case d = 2, so that ||T|| sep = 1 and we get, for the genuine four mode 
entanglement, 

£=||t||-1 (26) 

where 



\ i,j,k,l=l 



with 

t ijk i = Tr[p a,, <g> aj <g> a k ® a t ] = (ip\(Ji <g> (Jj <g> a k <g> ai\ip). (28) 

where {ai} i = 0,1,2,3 are the generators of the SU(2) group (Pauli operators). The 
resulting entanglement in \ip') is 



E g {\i//)) = ~ (-6 + ^88 + 64a 2 + 32/3 + 64/3 2 + 10a 2 f3 2 + /3 4 ) 

_4 (4/a - 2ga(l + /?) + /a/3(a 2 - /3 2 ) + 4ar/(l + /?)) e 
(-6 + ^88 + 64a 2 + 32/3 + 64/3 2 + 10a 2 /3 2 + /^ 



0[e 2 ] 

(29) 



where the first term gives the entanglement E(\ip)) for the state ip as defined in Eq. (l2~2l) . 
For this partition, the operations on a single mode are the only local operations, while all 
others are non-local. Therefore, the terms Thi and 777,3 are the only local interactions. 



12 



Therefore, we expect that the four mode genuine entanglement should not depend on T 
or 7 to the first order in e, which is the case, as seen from Eq.( l29l . 

Next, we consider the partition consisting of two subsets, each containing two modes 
on each site, {A \,A \.} and {B f, B |} (site partition). Thus we have two subsystems 
with d = 4 corresponding to a SU(A) <8> SU(4) qudit system. Further, n = 2 giving 
m = (2L/n) = 2 so that the geometric entanglement is 



E s (\^}) = \\r\\-Q (30) 



where 



\T\ 



15 



\ j,fc=i 



with 



Kjk = <^|Aj- <8) AfcIV) 

where A j ; j = 1, . . . , 15 are the generators of 5*f/(4). The entanglement of \ip') in Eg. (125 
is then given by 

E s (\i)')) = - (-18 + ^208 + 136a 2 + 9a 4 - 32/3 + 104/3 2 + 34a 2 /3 2 + 9/3 4 



(31) 



16(-/a + /a/3(a 2 -/3 2 -2))e + Q 2 



3(^208 + 136a 2 + 9a 4 - 32/3 + 104/3 2 + 34a 2 /3 2 + 9/3 4 ) 



According to the 'site partition', in addition to the operations on single modes, the 
operations on the pair of modes having the same site label are also local. Therefore, the 
resulting entanglement cannot change under the intra-site operations in the Hamiltonian, 
namely the q term, the 7] term and as before, T and 7 terms. Thus, to the first order in 
e, the entanglement is expected to depend only on the non-local part of the Hamiltonian, 
that is, on the / parameter. From Eq. 1 13 lft we see that this is the case. The inter-site 



14 



entanglement quantified using Von-Neumann entropy has been reported earlier J4|. 

Thus we see that the geometric measure has the capability to quantify the genuine 
multi-mode fermionic entanglement as against the mainly bipartite entanglement reported 

13 



in the literature. Also, the geometric entanglement measure, for the given partition of 
modes, shows the correct behavior under local and non-local unitary operations. 

B. Upper bound for the geometric measure in a four mode system. 

We discuss here some upper bounds that one can find for the entanglement for the 
four mode system and compare with existing results obtained from other measures. We 
find that the general state which leads to a maximum inter-site entanglement computed 
via the geometric measure also leads to a maximum for the von-Neumann entropy. 

We treat the entanglement for the given partition to be the function of the coefficients 
of the general state ]?/>) € #4(2), namely, 



E c *i*> ( 32 ) 



fe=3,5,6,9, 10,12 

where each ket is labeled by the (four bit) binary representation of k. We then maximize 
E(\ip)) with respect to the coefficients c^. 

For the site partition (see above) the entanglement (Eq. ( 1301) and two equations follow- 
ing it) is a function of the coefficients Ck in state \ip) G H^) given in Eq.(j32l. We find 
that the the maximum value of the entanglement is given by 

E max = 1.74593 (33) 

The Schmidt canonical form (with respect to the H^{2) basis) of \ip ma ,x) corresponding to 
this E mn:r . is 



■'max 



\ ( £ w) 

\/c=3,5,10,12 / 



\^ma X ) = -\ >. |*)| (34) 



or. 



\ ( E 1*)) 

\fc=3,6,9,12 / 



|Vw> = - I > , \k) ) (35) 



where each ket is labeled by the four bit binary representation of k. Also, both these forms 
lead to the von-Neumann entropy = 2 which is the maximum possible [131 ]. 
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For the four mode entanglement (Eg. (1261) ) we find that entanglement is maximum for 
the state given by 

c 3 = 0.06701 - 0.33751z c 5 = 0.16442 - 0.51281a c 6 = -0.1006 + 0.2854i 

c 9 = 0.29967 - 0.04206i c w = -0.53522 + 0.05955i c 12 = -0.34410 - 0.00118* 

(36) 

and the maximum value is found to be 

However, we do not have any entanglement measure to compare with the geometric mea- 
sure. We also note that, for four modes, no canonical form like the Schmidt or Acin 
canonical form (for two and three modes respectively) is available. 

V. APPLICATIONS 

We now consider some correlated fermionic lattice models and discuss multi-mode 
entanglement in these models using the geometric measure. 

A. Hubbard dimer. 



The Hubbard dimer model is a simple model for a number o 



ing the electrons in a H 2 molecule, double quantum dots, etc 
can be written as 



flE 



physical systems, includ- 



30j. The Hamiltonian 



H = -tJ2 ( C A* C B* + c ] Ba c A( /j + U J2 %t«jU ( 38 ) 

<T=t4 j=A,B 

where A, B are the site labels and \ and \, are spin labels, t is the hopping coefficient 
measuring hopping between two sites while conserving spin and U quantifies Coulomb 
interaction between fermions on the same site. By varying (^) we can vary the relative 
contributions of hopping and Coulomb mechanisms. 

15 



QEi 



The ground state of the system at zero temperature can be easily obtained as |4], |6| 

|V>o) = NG \vac) (39) 
where N = (V'olV'o) -1 ^ 2 is the normalization factor and 

G = c M c Ai + c m c Bi + a f — J [c M c Bi - c Ai c Bt J (40) 



with a(x) = x + yi + x 2 . By mapping to i?2L via Eq.(|2]) we get, 

\vac) -> |0) <g> |0) ® |0) ® |0) = |0000) (41) 

while mapping between the operators, (Eq.flHJ) and Eq.flSJ)) gives 

c At ^a) ®I 2 ®h®h c Al ->■ I 2 ®a) <g> i" 2 <g> I 2 
c B1 . ->■ J 2 <g> J 2 ® a f (g) J 2 c^ ->■ J 2 ® ^2 ® ^2 ® a f 



(42) 



The normalized ground state can be expressed in the qubit space as 



\ip ) = = {|1100) + |0011) + ot|1001) - a|0110)} (43) 

The four-partite entanglement in any state \ip) can be calculated by using 
Eqs.((T26]),(j27D,(lZHD) as, 

E b =\\t\\-\ (44) 

where 



:-! 



s E $■« ( 45 ) 

\ i,j,k,l=l 



with 

%w = (ipWi ® crj ® o-fc ® CT/IV'}- (46) 

The ground state entanglement can be then calculated to be 



B »=(TT^)V 1 + ^ 2 + a4 - 1 (47 > 
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We plot the four-partite entanglement as a function of U and t (Fig. 1(a)) and as a function 
of a (Fig. 1(b)). The entanglement is seen to monotonically increase as a function of a, 
saturating at large values of a to the maximum value 2. The saturation to the maximum 
value can be obtained either for very large values of U or very small values of t. We can 
interpret this result in the following way: since the total particle number is fixed to be 2, 
the four mode entanglement essentially measures the correlations between the spins. The 
entanglement increases as a function of a because the spin correlations increase with a. 

U 
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Fig. 1 Four-partite entanglement for the Hubbard dimer (at half filling) as a function of 
(a) (color online) U and t (both in energy units) and (b) as a function of a (a > 1). 



We can also calculate the bipartite entanglement between sites A and B using the 
geometric measure Eqs. (115|16|T71) considering the partitions to be {{A | A |}; {B f B \. 

}} 

E, = llrll - HrlL«, (48) 



I sep 



where 
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I A,: ® An 



(49) 



m/2 



and | \r\\sep = ( ^ ) • Here As are the generators of S77(4), there are m = 2 partitions 
(same as the number of sites) and each partition has dimension d — 4. This leads to an 
inter-site entanglement of the form 

2 



£, 



;i + « 2 ) 



V130 4 + 34a 2 + 13-6 



(50) 
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The bi-partite entanglement between sites A and B was calculated earlier using the von- 
Neumann entropy 4] 



E 



VN 



l + o 2 



log 2 [2(l + a 2 )] -a 2 log 2 



or 



2(1 + a 2 ) 



(51) 



We plot the inter-site entanglement (the von-Neumann entropy is also plotted for com- 
parison) as a function of a in Fig.2. It is seen that both measures show qualitatively 
similar behavior, i.e, a monotonically decreasing entanglement as a function of a saturat- 
ing at very large values of a. The entanglement between the sites A and B decreases as a 
function of a because with increasing on-site repulsion U, the four dimensional local state 
space at each site gets reduced to a two dimensional local state space |4| due to a sup- 
pression of charge fluctuations or in other words, as a — > oo the <S7/(4) eg) SU(4:) partition 
goes over to a SU(2)®SU(2) partition. We have explicitly checked that the entanglement 
obtained in the a — > oo limit matches with that obtained for the SU(2)®SU(2) partition. 





Fig.2 (a) The bi-partite entanglement between sites A and B calculated with the geometric 

measure as a function of a (a > 1) for the Hubbard dimer at half filling, (b) The 

corresponding von-Neumann entropy as a function of a (a > 1). 

We can also discuss bi-partite entanglement with unequal partitioning. Consider four 
modes partitioned into two subsets containing one and three modes respectively. This 
would correspond for example to one observer controlling a register which measures the 
occupancy of the spin up at site A and the other observer controlling a register which 
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measures the occupancy of the spin down at site A as well as the occupancy at site B. 
The bi-partite entanglement in this can be obtained as discussed in Sec.l lHIBl) . In the 
present case, the two partitions have dimension as d\ = 2 and d,2 = 8 respectively This 
gives, for the entanglement, 



E 



ii:<> 



\ j=\ k=\ 



or, using \ip) in Eq.( l4"3"|) . 



E = 1.6367 



(52) 



(53) 



Interestingly, the entanglement is independent of a. It turns out that this is the maximum 
value possible for the entanglement (Eq. (l52]) )( we checked this by maximizing the entan- 
glement given by Eq.f l52|) as a function of the coefficients in the general state Eq. (l32[) . We 
have also checked that von-Neumann entropy in this case is also independent of a and 
has the maximum possible value, i.e. 2. 



B. 



Three electrons on three sites. 



We next consider the Hubbard trimer, i.e, electrons on three sites with the sites A, B, C 
with periodic boundary conditions. Ignoring the chemical potential , the Hamiltonian is 



Htf) 



-t 



Yl ( c it c i+!t + c ]i c J+U + h - c ) ~ P Yl ™it%l 

J=A,B,C j=A,B,C 



(54) 



where t > is the hopping parameter and j3 



We begin our analysis as earlier by 



mapping from the fermionic to the six qubit space corresponding to (A^ A \. B \ B \, 
C f CI). In the qubit space, the basis respecting number super-selection rule is given 

by 

{\k)} k = 7, 11, 13, 14, 19, 21, 22, 25, 26, 28, 35, 37, 38, 41, 42, 44, 49, 50, 52, 56 (55) 

where the twenty basis states are labeled by the six bit binary representation of the k 
values. We numerically diagonalize the Hamiltonian H(f3) in Eq.( )54|) in the basis given by 
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Eq. (|55p at different values of as increases from zero. The total spin quantum number 
S and the z component of the total spin Sz commute with the Hamiltonian and can 
therefore be used as good quantum numbers to characterize the states. The ground state 
has total S value = 1/2. The triangular geometry of the three site model (with periodic 
boundary conditions) also leads to an additional symmetry under reflection about one 
of the medians of the triangle. This leads to a two fold degeneracy for the ground state 
(for a fixed Sz value). Since these symmetries are preserved even in the presence of the 
interaction U, the ground state remains two-fold degenerate for all values. 

The entanglement in any state can be calculated in a similar manner as shown previ- 
ously (Sec JVAj) . We now show the results of our calculations for the entanglement in one 



of the ground states. 

The six-mode entanglement as a function of the interaction parameter is shown in 
Fig.3a while the tripartite entanglement (between the sites A, B, C) is shown in Fig.3b as 
a function of 0. The tripartite entanglement between the sites A, B, C is seen to decrease 
with increasing - we interpret this result in a similar way as that for the dimer as due to 
the fact that the local state space at each site decreases with increasing 0. We find that the 
six-mode entanglement increases with 0, saturating at large values, however, the behavior 
is not monotonic. Such non- monotonic behavior as a function of is also shown by the 
bipartite entanglement between sites A and BC calculated using the geometric measure 
as well as the von-Neumann entropy (Figs. 3c and 3d). 

We can understand the non-monotonic behavior in the following way: the bi-partite 
entanglement measures the entanglement between the the sites A and the sites BC. With 
small increase in from zero, there is an increase in the spin fluctuations between B 
and C which shows up as an initial increase in the entanglement between A and BC. 
However, with further increase in 0, the charge fluctuations get completely suppressed 
leading to an asymptotic behavior similar to that for the dimer. One can understand 
the non-monotonic behavior of the six-mode entanglement in a similar way. For small 
0, there is a larger correlation between two sites which leads to a decrease in the overall 
entanglement, however with further increase in 0, the contribution to the entanglement 
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is solely due to spin fluctuations which increases with (3 leading to the observed increase 
in the entanglement as well. 
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Fig. 3 (a) Six-partite entanglement for the Hubbard trimer as a function of f3 at half filling, (b) 
Tripartite (site) entanglement as a function of j3 for the Hubbard trimer at half filling, (c) The 

bipartite entanglement between site A and sites BC (see text) as a function of f3 for the 
Hubbard trimer at half filling and (d) The corresponding von-Neumann entropy as a function 

of (3 for the Hubbard trimer at half filling. 



We have also calculated the upper bounds for the entanglements in various partitions in 



the same way as in Sec JIVl These are E ] 



max 
9 



4.42218, E% ax 



4.15105, E™ ax = 6.08767. 



VI. SUMMARY 

We have proposed a multipartite entanglement measure for N fermions distributed 
over 2L modes (single particle states). The measure is defined on the 2L qubit space 
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isomorphic to the Fock space for 2L single particle states. The entanglement measure is 
defined for a given partition of 2L modes containing m > 2 subsets, using the Euclidean 
norm of the m-partite correlation tensor in the Bloch representation of the corresponding 
multi-mode state, viewed as a m-partite state (see sec. IV and V). The Hilbert spaces 
associated with these subsets may have different dimensions. The quantum operations 
confined to a subset of a given partition are local operations. This way of defining en- 
tanglement and local operations gives us the flexibility to deal with entanglement and 
its dynamics in various physical situations governed by different Hamiltonians. We note 
that the concept of locality for indistinguishable fermions is distinct from that for distin- 
guishable particles. In the latter case, locality applies to spatially separated subsystems. 
However, spatially separated fermions become distinguishable. We have shown, using a 
representative case, that the geometric measure is invariant under local unitaries corre- 
sponding to a given partition. As an application, we have also considered some correlated 
electron systems and demonstrated the use of the multipartite measure in these systems. 
In particular, we have calculated the multipartite entanglement in the Hubbard dimer 
and trimer (at half filling). We find that the bipartite entanglement between sites as com- 
puted with the geometric measure has a qualitatively similar behavior as a function of the 
interaction as that of the conventional von- Neumann entropy. We have also calculated the 
four (six)- partite and the two (three) site entanglement for the Hubbard dimer (trimer). 
We find that the multi-partite entanglement gives complementary information to that of 
the site entanglement in both the cases. 

Although the entanglement measures given in this paper have been mainly applied 
to the study of the multipartite entanglement structure in Hubbard dimers and trimers, 
these measures are completely general and can be applied to other fermionic systems. We 



have also shown 
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20| that, viewed as a measure on qubit space, this measure has 
most of the properties required of a good entanglement measure, including monotonicity. 
To the best of our knowledge, this is the first measure of multipartite entanglement in 
fermionic systems going beyond the bipartite and even the tripartite case. Further, in 
this paper we have restricted to applications involving ground states of 2L mode systems 
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with 2L < 6. It is straightforward to extend the calculations for large number of modes 



(2L > 6) except for the length of the computation. We have shown earlier that 19|, |20] 
for antisymmetric states, the computational complexity of computing the measure goes 
polynomially with the number of parts of the system (here, the number of partitions of 
2L modes). It would be interesting to extend these calculations to larger system sizes 
where new and interesting results might be expected - we plan to do this in a separate 
study. 
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